Bond-versus-site doping models for off-chain-doped Haldane-gap system Y 2 Ba Ni O5 
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Using the density matrix renormalization-group technique, 
we calculate the impurity energy levels for two different ef- 
fective models of off-chain doping for quasi-one-dimensional 
Heisenberg chain compound Y2 Ba Ni O5: ferromagnetic 
bond doping and antiferromagnetic site spin- 1/2 doping. 
Thresholds of the impurity strength for the appearance of 
localized states are found for both models. However, the 
ground-state and low-energy excitations for weak impurity 
strength are different for these two models and the difference 
can be detected by experiments. 



In the last fifteen years many theoretical and exper- 
imental efforts have, been devoted to testing Haldane's 
original conjecture o that the excitation spectrum for in- 
teger spin Heisenberg antiferromagnetic (AF) chains has 
a finite energy gap above the ground state while it is gap- 
less for half-integer cases. The experiments on NENP and 
Y2 Ba Ni O5 have shown explicitly the gap's existence. 
The gap value can also be estimated by numerical stud- 
ies, e.g., for isotropic spin-1 AF Heisenberg chain, a gap 
A » 0.41J (where J is the exchange integral) has.be 
obtained by density matrix renormalization-group 
(DMRG) technique as welL-as by exact diagonalization 
with proper extrapolation o. However, the knowledge 
of these spin chains is still far from complete, especially 
regarding the impurity effects. Recently, the interest in 
the impurity problem has been revived due to the doping 
experiment performed on the linear chain compound Y2 
Ba Ni O5 0' O. Several theoretical studies have been 
carried out to study the doping effect tjQ. However, 
the off-chain doping (when Y 3+ is replaced by Ca 2+ or 
Mg 2+ ) has not been fully studied. This is the subject of 
the present paper. 

Y2 Ba Ni O5 contains one-dimensional chain structures 
composed of oxygen octahedra with Ni 2+ ions at the cen- 
ter. The inelastic-neutron-scattering (INS) experiments 
on single crystals of this charge transfer insulator show 
the chains are highly one dimensional withj-very weak in- 
terchain coupling J2 (J2/J < 5 x 10~ 4 ), V3 where J is 



the intrachain coupling, and the Haldane gap at q = tt 
is almost isotropic, so it is an ideal substance to con- 
front experiment. Furthermore, the doping of this com- 
pound can be in-chain as well as off-chain. The nonmag- 
netic Zn 2+ ions substituting Ni 2+ will effectively severe 
the chain and produce free s=l/2 spins at the two ends, 
which can. be explained by the valence bond solid (VBS) 
model. O Though the specific-heat measurements dis- 
playing Schottky anomaly H seemed to contradict this 
picture, a recent DMRG numerical study shows this con- 
flict can be resolved— jpy taking into account the chain 
length distribution. O On the other hand, the off-chain 
substitution of Y 3+ by Ca 2+ or Mg 2+ will induce holes 
along the chain. The INS experiments show the holes ate 
localized mainly on the oxygen sites within the chains, O 
between the nearest two Ni 2+ ions. Since the coupling of 
the two nearby Ni 2+ ions is mediated by this oxygen, the 
doping affect can be modeled in two different ways: one 
is that the coupling of the two Ni 2+ is modified (bond 
doping), the other is that an s=l/2 spin will be localized 
at the oxygen site, and the localized spin will be coupled 
to the two Ni 2+ (s=l/2 site doping). The Hamiltonian 
of these two effective doping models can be written as 
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respectively. Here the periodic boundary condition 
(PBC) is adopted, N is the chain length, S denotes the 
s=l spin; and in Eq.(^|) , cfq represents the s=l/2 impu- 
rity , and here we have assumed that the coupling of the 
two Ni 2+ ions around the impurity are destroyed com- 
pletely. 

From the site-doping point of view, the induced s=l/2 
spin couples antiferromagnetically with Ni 2+ on both 
sides of it. Pictorially, the spin direction of the two Ni 
spins around the impurity site will be antiparallel to the 
direction of the localized s=l/2 spin, that is, the two 
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s=l Ni spins tends to be in the same direction. So the 
effective coupling of the two site will be ferromagnetic 
instead of antiferromagnetic. Generically, in the Hamil- 
tonian p), the coupling J' will be negative, while for the 
site doping model in Hamiltonian (||), J" will be positive. 

Several numerical studies haffipb een . carried out for the 
two Hamiltonians (@) and (|), Mlitj but nearly all the 
calculations were done under the assumption that the im- 
purity bond remains antiferromagnetic. In this commu- 
nication, we use the DMRG algorithm to calculate the 
energy spectrum for the two doping models (Q) and (||), 
especially for ferromagnetic bond doping model. If J' 
(negative) is weak enough in Hamiltonian (jl|) , an in-gap 
bound state is localized around the impurity site. There 
is a threshold value J' cl above which there is no in-gap 
level or corresponding localized state. We also calculate 
the spectrum for the impurity state for antiferromagnetic 
J" in Hamiltonian Combined with previous studies 
of other authors ExJ , the two Hamiltonians with J' and 
J" AFM as well as FM coupling give a complete descrip- 
tion for the spin-1 Heisenberg chain doping models. Ex- 
perimental data may distinguish which of the two models 
is more appropriate in describing the off-chain doped Y 2 
Ba Ni O5 compound. 

The DMRG scheme we used is similar to the one used 
in Ref. |)|, we keep 250 states up to chain length 50 
throughout all the calculation, and the largest truncation 
error in our calculation is of. the order of 10~ 6 . At last 
the Shanks transformation D is used to obtain the final 
result. This transformation has been proved to.|be. valid 
for spin-1 Heisenberg chain by several authors ocm, and 
we can trust the accuracy of this scaling approach. Before 
presenting the numerical results, a physically intuitive 
analysis of the two doping models will be useful. 

First, for Hamiltonian (|l]) of the bond doping cases, re- 
sults in two limiting cases (J' — * 0~ and J' — > —00) can 
be obtained within the framework of the VBS picture. 
When J' = 0, the system becomes an open Heisenberg 
chain, and the ground state will be fourfold degenerate, 
where the two edge s = 5 spins couple into a singlet and 
a triplet. The lowest excited state with 5f ot = 2 is sep- 
arated from the ground state by the Haldane gap, and 
the continum spectrum starts there. When the impurity 
bond intensity increases gradually ( | J' | < < J, J' < 
in perturbation regime), the two edge spins will tend to 
form a triplet, while the singlet state will enter the Hal- 
dane gap and become an in-gap state between the triplet 
ground state and the continuum spectrum. The jyertur- 
bative result obtained by S0rensen and Affleck (SA) 
is valid also for the negative J' case. On the other hand, 
for the strong impurity doping strength, J' — > —00, the 
two s = 1 spins at the ends of the impurity bond cannot 
be reduced to s = | spins anymore. We must consider 
the reduced four sites model near the impurity bond and 
write down the effective Hamiltonian as 

H = ^e//o"i • Si + J'Si • S 2 + J e ffS 2 ■ er 2 j (3) 
where Si and S2 are s — 1 spins, and o\ and 02 are 



s = \ spins due to the VBS picture (the remaining part 
of the periodic chain can be viewed as a length L - 2 open 
chain) . 

In fact, the model Hamiltonian (||) can account for the 
essential features of the Hamiltonian (0) for all values of 
J', from positive to negative beyond the assumption that 
Jeff is positive. Its spectrum contains two singlets, four 
triplets, there quintets and one septet. Out of the ten 
states, only a few are important to our discussion: the 
ground state for J' > and that for J' < 0, the quintets, 
which are the low excitation states for some values of 
J' < 0. The energy of these states versus the impurity 
coupling J' (in units of Jeff) is shown in Fig. |l|. 
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FIG. 1. The low-energy spectrum of the reduced bond 
doping Hamiltonian (^). 

For the J' < case on the left part of Fig. |l|, we 
see that the ground state is a triplet for all J'. When 
J' — * -co, the first excited state is a quintet instead 
of a singlet. This is an interesting phenomenon for the 
ferromagnetic impurity doping case, that is, there will 
be a possible energy crossover between the singlet and 
quintet for larger \J'\. Since the energy difference of this 
quintet to the ground state is nonzero, we can expect 
either one bound state for small J e ff (more if J e // is 
small enough) or no bound states for larger J e // within 
the bulk gap A = 0.41048(2). We will show later that 
for large | J'\ there is no such crossover from our DMRG 
calculation, and there is no in-gap states for long chains. 
Therefore, J e ff is large enough. Since we know very 
little about J e // at the present stage, it is impossible 
to conclude which case will occur without more detailed 
computation. 

For J' > 0, the AFM doping case on the right part of 
Fig. [I], only two states are important: the singlet ground 
state and the first triplet excitation state. For J' = and 
J' — * co, the energy difference is zero, so the gap state 
can be expected for both weak and strong impurity bond 
strength. Detailed numericalrf esults for this case can be 
found in previous studies. EmS 
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FIG. 2. The energy spectrum of the bond doping Hamil- 
tonian (jjj) when the impurity bond is ferromagnetic. The dot- 
ted line indicates the Haldane gap of a pure spin-1 Heisenberg 
chain. 

Now we present DMRG results for Hamiltonian (|l|) 
with FM bond doping. The dependence of the relative 
energy of the lowest excited state with respect to the 
ground state on the impurity bond is shown in Fig. ||. 
From the figure, we see clearly that a threshold J' c ex- 
ists above which the in-gap state disappears. For weaker 
impurity bond, the impurity level increases linearly with 
the bond strength, which agrees with SA's perturbative 
argument for the weak coupling limit. The relation 
between the splitting of the ground state and the sin- 
glet excitations is AE = —a 2 J', where AE= 0.112J for 
J'=-0.1J, which is also in excellent agreement with the 
theoretical work for the AFM case 0. When the im- 
purity bond strength becomes bigger, the energy of the 
in-gap state deviates from linear and enters the contin- 
uum spectrum at last at a critical point for bond strength 
near J' = — 0.7 J. The simple four-site model based on 
the VBS picture discussed above shows a crossover be- 
tween the singlet state (in-gap state for some values of J' 
) and the quintet state, which remains at the bottom of 
the continuum spectrum, we did not find such a feature 
in the calculation, so there is no energy crossover taking 
place, this is due to the interaction between the singlet 
state and the continuum spectrum, which keeps the sin- 
glet state as the lowest-energy excitation state for longer 
chain length. 

The local bond energy (Si • Si+i) m p — (Si • Sj+i)i- for 
different impurity bond strengths is presented in Fig. || 
(here m represent the total spin and p is the parity of 
the state). The to p _— -2 + state corresponds to the one- 
magnon excitation V3x3 . It change very little when the 
impurity bond strength changes from very weak (-0.1J) 
up to very strong (4.0 J) (this is different from the posi- 
tive J'xase, because for J'— 1.0 J, this is a two-magnon 
states HE3). While the ra p = + state is an edge exci- 
tation localized around the impurity bond for weak J'. 
When J' becomes bigger, this excitation is delocalized 



and shows a single-magnon-like behavior in the strong 
coupling region. This is also consistent with the absence 
of the in-gap impurity states in this region. 




Fig.3-(a) : J =-0.1J 
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FIG. 3. The local bond energy {Si-Si+i) m p — (S i -S i+ i) 1 - 
for different ferromagnetic doping bond J', where m p equals 
0+ and 2+. 

We also calculate the ground state and the first ex- 
cited state (the subgap state) pair correlations of the site 
nearest to the impurity bond for different J' , that is, 
(Si ■ S50), which is presented in Fig. |]. For J' = 0, the 
open chain case, we see that for i > 25, the pair corre- 
lations for the two states are in phase, but for i < 25, 
they are out of phase. This phase difference which keeps 
the orthogonality of the two states is due to the differ- 
ent coupling of the two chain-end 1/2 spins. In contrast, 
for strong J 7 , the correlations are almost identical except 
those i in the middle of the-ftbain where the one-magnon 
amplitude is the biggest. Hi The correlations for the 
two states are both independent of J' for i > 25, while 
for 5 t z ot =l, the ground state, the correlations for i < 25 
will become larger, and for S* ot = 0, the first excited 
state, the correlation will change its phase when increas- 
ing J', the orthogonality of the two states will mainly 
come from those sites in the middle of the chain, this 
also means that the edge excitation (near the impurity 
bond) when J' is weaker will become the bulk excitation 
(far away from the impurity bond) when J' is stronger. 
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results are consistent with the previous calculation. 
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Fig.4-(b) : J =-0.4J 
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FIG. 4. The pair correlations (Si ■ S50) of the ferromag- 
netic bond doping for the ground state 1 _ and the first exci- 
tation state + . 

Secondly, we analyze the spectrum of Hamiltonian (||). 
The impurity coupling J" is positive. Let us reduce the 
Hamiltonian (0) to a three sites effective Hamiltonian 
around the impurity s = 1/2: 



H = Jeff(<ri • cr + (To ' Cjv). 



(4) 



Here ci and ctn are the effective s = 1/2 edge spins from 
the bulk of the chain, and <To is the impurity spin. As 
we analyzed above, here we have J e // > 0. In the per- 
turbation regime {Jeff ~ 0), we know that Hamiltonian 
(|[) has a doublet ground state, a doublet first excited 
state, and a quartet highest excited state. tJ For the orig- 
inal Hamiltonian (|^), as J e // increases, we could see that 
the doublet ground state will not change, while the two 
in-gap states given by the reduced Hamiltonian increase 
linearly with J e ff- By DMRG calculation, we find J" in 
Eq. (0) has two critical impurity coupling J" c and J^,, 
which correspond to the points where the doublet and 
the quartet enter the continuum spectrum, respectively. 

For the site doping Hamiltonian (||), we consider both 
the ferromagnetic coupling (J" < 0) and the antiferro- 
magnetic coupling case (J" > 0). The ground state is a 
quartet for the former while a doublet for the latter. Two 
in-gap states are found in the weak coupling regime in 
each of the two cases, and its energy spectrum is shown 
in Fig. [| The thresholds where the in-gap state en- 
ters the continuum are J" c = —0.25 J, J'{ c = —1.4 J and 
J" c = 0.3J, J% c = 0.5J for the ferromagnetic and an- 
tiferromagnetic impurity coupling, respectively. These 
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FIG. 5. The energy spectrum of the site doping Hamil- 
tonian (|^). The dotted line indicates the Haldane gap of the 
pure spin-1 Heisenberg chain. 

In summary, we have considered two kinds of effective 
models for off-chain doping in a one-dimensional Heisen- 
berg chain using the DMRG algorithm. Thresholds of 
the impurity bond strength are found for both ferromag- 
netic bond doping model and antiferromagnetic s = 1/2 
site doping model. The ground states of the two kinds 
of effective models are of total spin 1 and total spin 1/2, 
respectively. The in-gap states that appear for the weak 
impurity bond are also different for the two models. A 
singlet impurity state for ferromagnetic bond doping and 
a doublet along with a quartet for AF site s = 1/2 dop- 
ing. We believe these properties of the spectrum for dif- 
ferent effective doping models can be distinguished by ex- 
periments. The doped holes in the off-chain doping case 
can contribute to the transport properties of the Hal- 
dane chaM systems. This issue was discussed by other 
authors O and is not addressed here. 
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